Step Repulsion Mediates Wandering on a Si(OOl) Vicinal Face 
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With a Si(001) vicinal surface in mind, we study step wandering instability on a vicinal surface 
with an anisotropic surface diffusion whose orientation dependence alternates on each consecutive 
terrace. In a conserved system step wandering takes place with step-up adatom drift. Repulsive 
interaction between steps is found indispensable for the instability. Monte Carlo simulation with a 
strong repulsive step interaction confirms the result of linear stability analysis, and further shows 
that in-phase step wandering produces straight grooves. Grooves widen as their amplitudes increase 
in proportion to the square root of time. 
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On a vicinal surface of a crystal, steps undergo two 
types of dynamical instabilities: wandering and bunch- 
ing. Step wandering is the instability for step deforma- 
tion along the step, and step bunching is the one for 
the inter-step distanceEl. Both instabilities are caused by 
some asymmetry in the surface diffusion field. There are 
many effects which cause the asymmetry, and the drift 
flow of adsorbed atoms (adatoms) is one of them. On 
Si(lll) and Si(OOl) vicinal faces, a direct electric current 
induces drift of adatoms and the instabilities have been 
observed under its application^. 

Si(OOl) surface is reconstructed and forms rows of 
dimerized atoms arranged in a 2 x 1 unit cell (Figure [I]). 
On the reconstructed surface, the adatom surface diffu- 
sion is anisotropic such that it takes place more easily in 
parallel to the dimer rows than in perpendicular. On a 
vicinal face terraces of different heights are bounded by 
steps. On consecutive terraces the dimer orientation al- 
ternates, and we call the 1x2 terrace and the 2x1 
terrace Tg. 

Due to the alternation of the orientation of fast surface 
diffusion on different terraces, conditions of the step in- 
stabilities for a Si(OOl) vicinal face differ from those for 
a Si(lll) vicinal face. Experimentally, bunching is ob- 
served on a (001) vicinal .face with a finite current irre- 
spective of its directipnErffl, and the step wandering with 
the step-up currcntu. .-Since the drift is believed to be 
parallel to the currentQ, the drift direction to cause the 
step wandering is opposite to that on a Si(lll) vicinal 
face. 

Theoretically, step bunching on a Si(001) vipinaLface 
is studied by a onc-dimcnsio>aal step flow modcB'BE£l and 
by Monte Carlo simulations^!]. When the alternation of 
anisotropic surface diffusion is taken into account, the 
step bunching instability is found irrespective.jOf.the drift 
direction, in agreement with the experimentstjQtj. On the 
contrary, there is no theoretical study on the step wan- 
dering on a Si(OOl) vicinal face so far, which we undertake 
in this paper. 

Atoms detached from steps migrate on terraces and 
attach to some steps. For the wandering instability, the 



adatom drift is necessary as well as the diffusive mo- 
tion. Evaporation and impingement are omitted. Steps 
are running parallel to the x-direction on average, and 
the positive y-direction is chosen in the step-down di- 
rection. The drift is assumed in the y-direction. With 
the anisotropy of the diffusion coefficient, the diffusion 
equation of adatom density c(x,y,i) is expressed as 

d t C = D X d 2 x C+Dydlc- fDydyC, (l) 

where dt, d x and d y represent the partial derivatives 
with respect to the subscripts, D x is the diffusion co- 
efficient in the x-direction, D y in the y-direction, and 
/ (= F/k-gT) the force to induce the drift devided by 
the temperature. The meaning of D x and D y depends 
on which terrace we are discussing, on T^ or Tb, and 
we come to this point later. For simplicity, we assume 
that the step kinetics is fast enough that the adatom 
density attains its local equilibrium value at each step: 
c| ± = c° eq [l + 0,0k + d y ()/k B T}. Here +(-) indicates 
the lower (upper) side of the step, Cg q is the equilibrium 
adatom density of a free straight step, f2 is the atomic 
area, /3 the step stiffness, n the curvature of the step, and 
C the step-step interaction potential. We assume that £ 
is a function of the step distance / in the y-direction as 
C = A(l+ 2 + IZ 2 ) with a positive constant A, correspond- 
ing to step repulsion. There is a more detailed model 
where £ is expressed by an integration of the force dipole 
along the steplj. But in the linear analysis the complica- 
tion is shown to be incorporated into the renormalization 
of the stiffness. Therefore, we use the simple form for ( 
here. 

By solving eq. ([!]) with the boundary conditions in the 
quasi-static approximation (d t c = 0), the adatom density 
c and then the adatom current j are determined. The 
step velocity is given by V — Vln ■ — > where n 
is the normal vector in the step-down direction. 

Due to the different orientation of the dimer rows on 
Ta and Tb, a set of diffusion coefficients (D x , D y ) corre- 
sponds different combinations; (D x ,D y ) = (D\\,D±) on 
Ta and (D x ,D y ) = (D±,D») on Tb, and D± is larger 
than Dm. Since the step Sb is rougher than S,4 on a 
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Si(OOl) vicinal face, step parameters arc different in gen- 
eral for the two types of steps, S^ and Sb, but for sim- 
plicity, we neglect these differences. 




FIG. 1: A Si(OOl) vicinal face. Short lines on terraces repre- 
sent dimers. 

On a flat vicinal face where parallel steps are arranged 
cquidistantly, the adatom concentration is homogeneous 
as Cg q . The step velocities Va and Vb of the steps Sa 
and Sb are calculated as Va = —Vb = —n.f(D^ — D\\)c® q . 
Since the step velocities Va and Vb are nonvanishing and 
opposite for a finite drift (F ^ 0), the flat vicinal face 
is unstable against step pairing. Without the repulsive 
step-step interaction, the adatom concentration is homo- 
geneous as Cg q irrespective of the step distance. The steps 
move with the velocities Va and Vb given above, and by 
coalescence the surface is covered by one type of terrace; 
for instance, T^ for positive F since D± > Dn. 

If the repulsive interaction is granted, the difference of 
terrace widths I a and Ib of terraces T^ and Tb causes the 
difference of the equilibrium adatom densities ca and c-q 
at steps Sa and S b ■ Then the steady state with vanishing 
step velocities can be established even with the drift as 

£>±(c A e^ -cb) = D||( CbC /'b _ CA ) 

e flA _ 1 e fl B - I ■ { > 

Since the average terrace width of the vicinal face is 
fixed to Z, wide and narrow terraces appear alternately 
as /a = I + 51/2 and Ib = I — 51/2, and the equilibrium 
densities also alternate as ca = Cg q — Ac/2 and cb = 
c" q + Ac/2, respectively. For a small drift fl and the 
strong repulsive interaction, the density difference Ac = 
cb — ca is approximated by Ac/Cg q w <rfl with a = 
(£)_!_ — D\\)/(D± + Dn), and the terrace width difference 
is 51/1 f» ksTafl 2 /12f2A Without repulsive interaction 
A = 0, the terrace width cannot remain finite. Also the 
above equation indicates that the deviation of the terrace 
width 51/1 is small under a very strong step repulsion. 

We now study the stability of this steady state with al- 
ternating terrace widths under the adatom drift. When 
the step repulsion is strong enough, the perturbations to 
step positions 5yA{x,t) and 5yB{x,t) of a wavenumber q 
can be shown to merge into the in-phase fluctuation of 
the same amplitude 5yA{x,t) = 5yB(x,t) = 5y q e u ' qt e ,J ' !x . 
With this deformation the terrace widths in the y di- 
rection remain Ia and Ib as before. For small q, the 



amplification rate is calculated to be 

LU q = a 2 q 2 - a^q 4 , (3) 

where a 2 = -fl(D±_ - Dp) Ac/2 and a 4 = 0, 2 c° q (D\\l A + 
D jJb) $ l '2/cbT '. Here we have assumed that the terrace 
widths are small as ql a,b *C 1. 

When the drift is in the step-down direction(F > 0), 
I A is larger than Ib, and cb is larger than ca- Then the 
quadratic term in q is negative and the steady state is sta- 
ble. When the drift is in the step- up direction(F < 0), 
ca is larger than cb, and the quadratic term is positive. 
Then the steady state with straight steps is unstable and 
in-phase step wandering occurs. Thus the first term in- 
dicates an instability of the steady state by changing its 
sign with the drift direction. Note that the wandering is 
caused by the difference Ac of equilibrium concentrations 
at the steps originated from the step repulsion. Thus 
the step repulsion is indispensable for the establishment 
of wandering instability on a Si(001) vicinal. Also note 
that there is no critical value for the step- up drift. If 
the drift changes from step-down to step-up, steps im- 
mediately shows the wandering instability, because the 
Gibbs-Thomson effect gives only a quartic term, higher 
than the destabilizing quadratic term. 

For a strong repulsive interaction, the difference of 
the step distance 51 is small and the wavelength of the 
most unstable mode is given by A max = 2ix 2a^/ 012 — 

2-ko~ x \j 2^l^/\F\. The characteristic wavelength is in- 
versely proportional to the square root of the external 
field E which causes the drift: A max oc l/yE. 

After the instability sets in, the step deformation am- 
plifies, and a nonlinear analysis is wanted. Assuming an 
in-phase motion of steps due to the strong step repul- 
sion, a heuristic argument is possible on the nonlinear 
evolution of the step wandering 5y(x,t) = r)(x,t). This 
simultaneously reveals the physical origin of the wander- 
ing. Since the number of atoms is conserved, the in-phase 
step motion is controlled by the adatom current in the x 
direction. The drift current has only the y component so 
that only the diffusion current determines the evolution 
of step fluctuation. 

There are two contributions of diffusion current. One 
is the current across the terrace induced by the differ- 
ence of the equilibrium densities at both ends of the 
terrace. When two bounding steps are tilted from the 
y direction with an angle 9 = t&n~ 1 (d x r]), the gradi- 
ent of the adatom density in a terrace Ta is Vc^ = 
(Ac/Ia cos#)(— s'mO, cos 8). The diffusion current is ob- 
tained as j A — (Du Act&n9/lA, — D±Ac/Ia), and the 
x component of the total flux on T^ is J x = j^Ia = 
Z)i|Actan# (one obtains the same result by solving 
(0)). Similarly the x component of the flux on Tb 
is J x = —D±Act&nd, and the average flux per step 

is Jl 1 ' = (J* + jZ)/2 = -(D x - D {] )Acd x r]/2. A 
steady state condition similar to eq. (||) with tilted steps 
gives the difference of the equilibrium adatom density as 
Ac = cr/ZCg q cos 2 6, and the extra factor cos 2 9 appears. 
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FIG. 2: Snapshots of the step wandering (a) with step-down drift at t « 5.0 x 10 , (b) with step- up drift at t « 2.5 x 10 5 , and 
(c) with step- up drift at t m 12 A x 10°. The number of steps is 32 and the system size is 256 x 256. 



When Ac < (i.e. F < 0), the flux is an increasing 
function of the slope and an instability is expected. 

The other is the diffusion current along the steps due 
to the chemical potential change with curvature. Consid- 
ering that the normal distance between steps for a tilted 
part is Ia,b cos 9, the flux (per step) passing through the 
two terraces induced by the chemical potential gradient is 
Jx 2) = -{lACOs9D ll +l B cos6D ± )c cq cos6d x n/2k B T. The 
chemical potential [i is determined solely by the Gibbs- 
Thomson effect Q,(3k and independent of the step repul- 
sion in the present choise of interaction £, since the ter- 
race width in the y direction is constant for the in-phasc 
modulation. More general expression given by Paulin 
and co-workerst3 might modify the following result quan- 
titatively but not qualitatively. 

With two contributions together, mass conservation 
leads to the following nonlinear time evolution of the in- 
phase step deformation 77(2, i); 

dtV = -^[^(^il + ^fff- 3 / 2 ^))], (4) 

where and g — l + (d x r]) 2 = sec 2 9. ct2 and 04 are those in 
eq. (Q). With 77 = 8y q e lqx+u}t the linear amplification rate 
ujq is recovered from cq.(||). Interestingly cq.([|) is the 
same as the npa jinear equ ation obtained for other con- 
served systemal3't3E3llall3 although the mechanism looks 
very different. 

The above theoretical analysis is now compared with 
the results of Monte Carlo simulations. For the simu- 
lation algorithm, consult Refs. |lll , p"7||l8| . We study steps 
with an average distance I = 8 on a square lattice sys- 
tem of size 256 x256 or 512 x 128. Length hereafter is 
measured in the unit of lattice constant a, and time in 
the unit of a 2 /D±. The parameters are so chosen as the 
equilibrium adatom density c° q = 0.18, the step stiffness 

/3/k B T = 0.13, D\\ = 0.5 and D±_ = 1.0. There is no 
extra energy barrier for the over-step diffusion. Kinetic 
coefficient is large enough so that the local equilibrium 
condition is valid. The strength of repulsive potential 



A/k B T — 46 is large enough to suppress step bunching 
in the following simulations. 

In Figure ^| we show snapshots of the step wandering 
under various drift forces; / = 0.1 for Fig. ||(a) and -0.1 
for Fig. |2|(b) and (c). Solid lines represent the step S,4 
and dotted lines represent the step 8b- We start the sim- 
ulation with an equidistant train of straight steps. With 
step-down drift in Fig. |^(a) (upward drift in the figure), 
steps remain straight. With step-up drift in Fig. g(b) and 
(c) step wandering occurs, in agreement with the linear 
stability analysis. Because the wandering is in phase, 
grooves appear parallel to the y-axis. 

Since \fl\ = 0.8 is not very small, we have to use a 
general formula for the wavelength of the most unstable 
mode, and it is obtained as A max sa 77, in good agreement 
with the result A « 64 in Fig. || (b). There, the wandering 
amplitude (the average step fluctuation width) is w 
14.5. In a late stage shown in Figure ||(c) the amplitude 
increases up to w « 37.2 when the wavelength of the 
grooves is about 85. Thus, the structure coarsens parallel 
as well as perpendicular to the steps. 

Recently, Paulin and co-workcrsB studied the step 
wandering induced by the Ehrlich-Schwoebel effect in a 
conserved system. They found perpetual enhancement 
of the wandering amplitude as w ~ t 13 with (3 ~ 1/2, 
irrespective of the step repulsion, but the coarsening of 
the wavelength of grooves took place only with a step 
repulsion. Although we take account of the long-ranged 
step repulsion only in the y-direction but not in the x- 
direction as Paulin et al. did, the entropic repulsion may 
have caused an effective repulsion in the x-direction and 
eventually the coarsening in our case, too. 

In Figure || open circles represent the time evolution of 
the wandering amplitude w for / = —0.1. The result is 
obtained by averaging over 10 runs of the size 512 x 128 
with 16 steps. In an early stage (t < 2 x 10 5 ), the step 
width increases rapidly. Then, the width enhancement 
slows down to i« ~ t 1 / 2 . The exponent is the same with 
the values obtained for the step wandering in other con- 
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FIG. 3: Time evolution of the step width (0)> w ~ t 1 ^ 2 , and 
the groove wavelength (□), A ~ t 1 ^ 6 . 

served systems0[liB00. The slowing down of the fluc- 
tuation amplification is attributed to the suppression of 
the diffusion current due to the narrowing of the terrace 
widtlM The groove wavelength A is obtained by count- 
ing the number of grooves for 10 samples, and is depicted 
by open squares in Fig. |^. The slow increase as A ~ t a 
with a = 0.17 ± 0.04 is consistent with the jape found 
with the use of a generalized version of eq. (4)E£l. 

In this paper, we studied the drift-induced step wan- 
dering on a vicinal face with an anisotropic surface dif- 
fusion whose orientation dependence alternates on con- 
secutive terraces. The step-step interaction is shown to 
play an essential role for the step wandering, since it not 



only prevents steps from coalescing but also creates the 
difference Ac of the equilibrium adatom density. The im- 
balance of the diffusion current between different steps 
induced by Ac causes the step wandering. Thus the 
step repulsion mediates steps to establish the asymme- 
try in the diffusion field. Due to this asymmetry, the 
step wandering occurs with step-up drift. The in-phase 
wandering leads to the formation of straight grooves on 
a vicinal-Jape, in accordance with other conserved sys- 
temsBBM. 

Recently, in an experiment by Nielsen and coworkers! 
a dimpled specimen was used and the step wandering on 
Si(001) vicinal face was observed with the application of 
a direct electric current. Near the bottom of the dimple, 
where the inclination is very small, the step bunching oc- 
curs irrespective of the current direction^. Fluctuation 
of bunches with the step-up current is larger than that 
with the step-down current. On increasing the inclina- 
tion, which means increasing the repulsive interaction, 
the in-phase step wandering occurs and straight grooves 
parallel to the current appear with the step- up current. 
The step wandering was observed in a range of inclina- 
tion angle between 0.08° and 0.5°. In Si(001) vicinal face, 
the direction of the drift of adatomS[4s believed to be the 
same as that of the electric currentQ, and our results of 
step wandering under a step-up drift qualitatively agree 
with the experiment!!]. 
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